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Single photon spectra in heavy-ion collisions at SPS energies are studied in the relativistic trans¬ 
port model that incorporates self-consistently the change of hadron masses in dense matter. We 
separate the total photon spectrum into ‘background’ arising from the radiative decays of tt® and 
r} mesons, and the ‘themal’ one from other sources. For the latter we include contributions from 
radiative decays of p, co, p', and ai, radiative decays of baryon resonances, as well as two-body 
processes such as tttt ^ py and vrp ^ Try. It is found that more than 95% of all photons come 
from the decays of tt® and p mesons, while the thermal photons account for less than 5% of the 
total photon yield. The thermal photon spectra in our calculations with either free or in-medium 
meson masses do not exceed the upper bound set by the experimental measurment of the WA80 
Collaboration. 
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I. INTRODUCTION 


One of the primary motivations of ultra-relativistic heavy-ion collisons is to create and study in the laboratory the 
properties of quark gluon plasma (QGP) as predicted by quantum chromodynamics (QGD). Experiments have been 
carried out at BNL-AGS and CERN-SPS energies, and will be continued at higher energies of BNL-RHIG and GERN- 
LHG colliders Q . The measurement of electromagnetic observables such as photon and dilepton spectra constitutes a 
major part of these efforts. The main reason for this is that the photons and dileptons do not suffer strong final-state 
interactions as hadrons do. They can thus be considered as ‘penetrating probes’ of the initial hot and dense stages 
where the QGP might be formed in these ultra-relativistic heavy-ion collisions ii- 

So far, the measurement of eletromagnetic observables has chiefly been carried out at CERN-SPS energies. Recent 
observation of the enhancement of low-mass dileptons in central heavy-ion collisions by the CERES and Helios-3 

collaborations ||^ has generated a great deal of interest in the heavy-ion community. Different dynamical models, such 
as hydrodynamical and transport models, have been used to investigate this phenomenon |1^-|^ . Calculations based 
on ‘conventional sources’ such as Dalitz decay and direct vector meson decay that account for the dilepton spectra 
in proton-induced reactions fail to explain the observed enhancement in heavy-ion collisions. Various medium effects, 
such as the dropping vector meson masses [p^[- |T^J^ , p5[ | as first proposed by Brown and Rho , the modification of 
rho meson spectral function and the enhanced production of p and/or p' [p^ , p0[ , have been proposed to 

explain this enhancement. 

Another piece of experimental data from CERN-SPS that have been discussed extensively is the single photon 
spectra in central S-|-Au collisions from the WA80 collaboration js^. In this experiment, 7r° and p spectra were 
measured simultaneously, so that their contributions to single photon spectra can be subtracted. It was found that 
the direct photon excess over these background sources is about 5% for central S-|-Au collisions pl|-^. Similar 
measurement of inclusive photon spectra has been carried out by the CERES collaboration . The results are very 
much in agreement with those of WA80, namely, the inclusive photon spectra can be basically explained by hadronic 
decay, particularly radiative decay of 7r° and p mesons. 

Many hydrodynamical calculations have been carried out to study photon spectra in these reactions ||35|-^. In 
all these calculations, it is found that the contributions from the QGP, if it were formed at all in these reactions, 
to the single photon yield is negligibly small. The major contributions to these so-called ‘thermal’ photons, namely, 
excessive photons over the background of 7r° and p radiative decays, come from hadronic interactions. 

In a number of papers |3qJ^,M, the absence of significant thermal photons has been interpreted as an evidence 
for the formation of quark gluon plasma. Without phase transition, the initial temperature of the hadronic gas was 
found to be about 400 MeV . This led to a large amount of thermal photons from hadronic interaction which 

was not observed experimentally. Including the phase transition, the initial temperature can be lowered to about 200 
MeV because of increased degrees of freedom. This lower initial temperature natually reduces the thermal 

photon yield to be in agreement with the WA80 data. There is, however, a major caveat in this type of analysis 
and reasoning, namely, the initial temperature in the hadronic scenarios depends sensitively on how many hadron 
resonances are included in the analysis. Indeed, in Ref. ||^, it was found that including a reasonable amount of 
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hadron resonances the initial temperature can be reduced to about 250 MeV. Furthermore it was shown in Ref. @] 
that, if all the hadron resonances with masses below 2.5 GeV are included, the initial temperature can be lowered to 
about 200 MeV. In both papers, it was shown that the WA80 photon data can be explained without invoking the 
formation of quark gluon plasma. 

There are, up till now, no detailed transport model calculations for photon spectra in heavy-ion collisions at CERN- 
SPS energies. In Ref. [^, results from the Hadron-String dynamics (HSD), that includes photons from the decay 
of 77 ', uj and ai, were compared with the data and were found to be far below the upper bound set by the WA80 
collaboration. However, photons can also be produced from the decay of rho mesons, baryon resonances, and two-body 
scattering processes such as tttt ^ py and 7rp ^ Try, which were neglected in Ref. Q. Whether the inclusion of these 
additional sources will bring the theoretical results close or even beyond the upper bound of the WA80 data needs to 
be studied. 

One of the purposes of this paper is thus to calculate the single photon spectra in central S-|-Au collisions based on 
the relativistic transport model and to compare with hydrodynamical results. The other purpose is to see whether 
the same model that explains the enhancement of low-mass dileptons, as shown in Refs. EHil, can also explain the 
lack of signal in the direct photon measurement. Explaining simultaneously these two correlated facts constitutes an 
important test for the consistency of the model. 

This paper is arranged as follows: In Section 2 we recall briefly the main ingredients of the relativistic transport 
model. In Section 3 we discuss various source for photon production, intending meson and baryon decay, as well 
as two-body scattering. We shall discuss radiative decay widths and photon production cross sections in two-body 
processes. The results are discussions are presented in Section 4. The paper ends with a brief summary in section 5. 


II. THE RELATIVISTIC TRANSPORT MODEL AND HADRONIC OBSERVABLES 


In studying medium effects in heavy-ion collisions, the relativistic transport model [|I3[-tI5| based on the Walecka- 
type model [|46| has been quite useful, as it provides a thermodynamically consistent description of the medium effects 
through the scalar and vector fields. In heavy-ion collisions at CERN-SPS energies, many hadrons are produced 
in the initial nucleon-nucleon interactions. This is usually modeled by the fragmentation of strings, which are the 
chromoelectric flux-tubes excited from the interacting quarks. One successful model for taking into account this 
nonequilibrium dynamics is the RQMD model j^ ] . To extend the relativistic transport model to heavy-ion collisions 
at these energies, we have used as initial conditions the hadron abundance and distributions obtained from the string 
fragmentation in RQMD. 

Further interactions and decays of these ‘primary’ hadrons are then taken into account as in usual relativistic 
transport model. We include non-strange baryons with masses below 1.72 GeV, as well as A, E, A(1405) and E(1385). 
For mesons we include tt, tj, p, uj, 77 ', oi, and </), as well as K and Ar*(892). Baryons are propagated in their mean 
fields, which are assumed to be the same for all non-strange baryons. The mean fields for hyperons are assumed to 
be 2/3 of that for non-strange baryons, based on the simple quark counting rule. Medium effects on pions and phi 
mesons are neglected, while those on kaons are always included, with their mean fields taken to be 1/3 of these for 
nucleons. For other non-strange mesons, we consider two scenarios as in Refs. Kl, namely, one with medium 
effects and one without medium effects. In our model, these mesons are treated in the constituent quark model. They 
thus feel only the scalar potential that shifts their masses, since the vector potential on the (constiutent) quark is 
canceled by that on the (constituent) antiquark. 

In addition to propagation in mean fields, hadrons also under stochastic two-body collisions. For baryon-baryon 
interactions, we include both elastic and inelastic scattering for nucleons, A(1232), V(1440) and fV(1535). Their 
cross sections are either taken from Refs. |48 4^ or obtained using the detailed balance procedure |^. The meson- 
baryon interactions are modeled by baryon resonace formation and decay. For example, the interaction of a pion 
with a nucleon proceeds through the formation and decay of A(1232), A(1440), ..., N(1720). The formation cross 
sections are taken to be of the relativistic Breit-Wigner form. The meson-meson interactions are either formulated by 
the resonance formation and decay when the intermediate meson is explicitly included in our model or treated as a 
direct elastic scattering with a cross section estimated from various theoretical models. Of particular importance for 
photon production are the pion-pion and pion-rho collisions, which are dominated by rho and oi meson formation, 
respectively. 

As mentioned in the Introduction, the observed enhancement of low-mass dileptons in heavy-ion collisions can be 
explained by including dropping vector meson masses. It is useful to see whether the inclusion of the dropping vector 
meson mass is also consistent with the real photon data from the WA80 collaboration. For this purpose we will also 
calculate photon spectra in the scenario of dropping meson masses. Following Ref. [^2|-Q, we extend the Walecka 
model from the coupling of nucleons to scalar and vector fields to the coupling of light quarks to these fields, using 
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the ideas of the meson-quark coupling model For a system of nucleons, pseudoscalar mesons, vector mesons, and 
axial-vector mesons at temperature T and baryon density ps, the scalar field (cr) is determined self-consistently from 
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where we have used the constituent quark model relations for the nucleon and vector meson masses i.e., m’^ = 
mjsf — gaW), m*p^ « mp^^ — {2/3)gcr{a), the quark structure of the rj meson in free space which leads to m* ~ 
mri — 0.45ga{a}, and the Weinberg sum rule relation between the rho-meson and ai meson masses, i.e, m^^ « 
mai — {2\^/3)ga{a). Recently we found that the use of a refined model, the effective chiral Lagrangian of [^, leads 
to essentially the same results for dilepton spectra Q]. 

This model describes the hadron observables quite well. In Fig. 1 we show the pseudorapidity distribtion of charged 
particles (upper window) and transverse momentum spectra of pions (lower window) for central S-|-Au collisions. They 
are found to be in good agreement with the experimental data from the WA80 and the CERES collaboration 

Comparisons with other hadron observables, such as proton rapidity distribution and transverse mass spectra 
can be found in Ref. |Q. 


III. PHOTON PRODUCTION: RADIATIVE WIDTHS AND CROSS SECTIONS 

The majority of single photons come from the decay of tt^ and g mesons at freeze out. The background photon 
spectra from their decay can thus be evaluated at the end of the transport model simulation by using the following 
branching ratio @ 




■0—►^--v — 


FO 


El = 0.988, R 




= ^EEUI ^ 0.3925. 


( 2 ) 


For the ‘thermal’ photon spectra we include the decay of p, oj, g\ and oi mesons, the decay of all the baryon 
resonances in our model, as well as two-body processes such as tttt —> py and np —> Try. The decay widths for 
p3 ^ TT+TT^y and p^ —> Tr^rr^y are taken from j^, which explains the measured width for p^ —> 7r+7r“y. The mass 
dependence of these widths are shown in Fig. 2. The open circle gives the experimental data for —> Tr+rr^y [ ^ . 
Note that in transport models, rho mesons can have different masses because of its broad mass distribution, so mass 
dependent widths are used. 

The omega meson radiative decay width was studied in Ref. ||^ ] based on the chiral Lagrangian that 

includes the Wess-Zumino anomalous term. The width was found to be proportional to the pion momentum in the to 
rest frame, namely. 


= C|p,|3 « 0.01316 |p^|3, (3) 

where |p 7 r| = (M^ — to^)/( 2M) and the coefficient C is determined from the measured width [^ . 

There exists only one experimental measurement of the ai radiative decay width, using the reverse process of 
oi production in pion-nucleus collisions |^l|. The width is found to be = 0.64 ± 0.246, which was used in 

Ref. in calculating ai contribution to photon production. In Ref. |^, an effective Lagrangian was proposed to 
calculate this width, which turns out to be about 1.4 MeV using the vector dominance assumption. This is close 
to the prediction of non-relativistic quark model of about 1.0-1.6 MeV [^. In view of large uncertainties, we use 
= 1 MeV as in Refs. |p^ , |T^ . According to Ref. [^, in the limit that pion momentum is much greater 
than its mass, the ai radiative decay width is given by 

= C\p^\^ « 0.01214|p,|5, (4) 

where the coefficient C is determined by using Ta^^Tr-fimai) = 1 MeV. 

Photons can also be produced from the decay of baryon resonaces. These contributions are usually neglected in 
hydrodynamical calculations ||^,^l|]. The radiative decay widths for baryon resonances included in our transport 
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model have been measured experimentally through photon-nucleon and/or pion-nucleon interactions j^ ]. We use 
these widths in our calculation of their contributions to single photon spectra. In addition, decay dominantly 
into Ay. But since the mass difference between and A is not very large, photons from E*^ decay usually have low 
transverse momenta, as will be shown below. 

For contributions from two-body processes, we include tttt —> py and irp — > Try, which are the most important ones 
in the temperature region relevant for SPS energies |^,^. These cross sections were first evaluated in Ref. with 
an effective Lagragian including tt, p, p, and lo mesons, which will be used in the present work. The cross section for 
TTTT ^ py is shown in Fig. 4, for three different rho meson masses. As the rho meson mass is reduced, the threshold 
and also the magnitude of the cross section are reduced. The cross section for np —>■ Try is shown in Fig. 5 for three 
different rho meson mass. When the rho meson mass is reduced, the available center-of-mass energy is reduced and 
the cross section increases. 

The importance of the oi meson was investigated in Refs. Since in our transport model the processes 

Trp ^ oi and oi — > Try are explictly included, we need not to include the effects of oi on the direct Trp —s- Try process. 
Otherwise there will be double counting. The effects of the oi on tt-tt ^ py were found to be appreciable for photons 
with energies greater than about 0.5 GeV |6|]. The contributions of tttt ^ py to these ‘hard’ photons are, however, 
extremely small (about one to two orders of magnitude smaller than those from Trp —> Try). This justifies our use of 
the cross sections given in Ref. that did not include the oi meson in the model. 


IV. PHOTON PRODUCTION: RESULTS AND DISCUSSIONS 

We present in this section the single photon spectra in central S-|-Au collisions at 200 AGeV. The acceptance of the 
WA80 collaboration is taken into account, namely, we include photons with pseudo-rapidity of 2.1 < rj < 2.9. In Fig. 
6 we show the background photon spectra from the decay of tt*^ (dotted line) and p meson (dashed line). It is seen 
that contribution from tt^ decay far dominates the background. 

The contributions to the so-called ‘thermal’ photons from the decay of mesons and baryons, and from two-body 
scattering are shown in Fig. 7. The uj radiative decay is found to be the most important source for photons with 
transverse momenta above about 0.5 GeV. This is in agreement with the finding of Ref. although our results are 
somewhat larger than those of Ref. This difference might be traced back to the difference in the ‘primary’ lo 

meson abundances in the RQMD and the HSD models. The contribution from the ai radiative decay is comparable to 
that from direct Tip —> Try. This is in agreement with the conclusion of Ref. |]^ that including ai the Trp contribution 
increases by about a factor 2-3 in the relevant temperature and photon energy region. 

The contributions from rj' and E° radiative decays are restricted to photons with transverse momenta below 0.5 
GeV, because the mass differences between hadrons in the intial and final states in these processes are small. Photons 
with transverse momenta below 0.2 GeV come chiefly from the decay of E° and tttt scattering. The reaction tttt — > py 
is endothermic, with most of the available energy going into the rho meson mass. This cross section actually diverges 
as the photon energy goes to zero. We have included a low-energy cut-off as in Ref. Q|. The choice of this cut-off 
parameter affects basically only the photons with transverse momenta below about 0.1 GeV, where no experimental 
data are available. 

In Fig. 8 we compare our results with the upper bound of the WA80 collaboration. More than 80% of thermal 
photons come from meson decay, of which more than half come from the w radiative decay. Overall, our results are 
well below the WA80 upper bound for photons with transverse momenta below 1 GeV. For higher tranverse momenta, 
our results touch the upper bound of the experimental data. Our results are actually very close to those of Ref. |^] 
with the parameter set named EOS H, that includes about the same number of hadron degrees of freedom as in our 
model. The conclusion is thus that the WA80 single photon data do not necessary imply the formation quark gluon 
plasma. Of course, the possibility of quark gluon plasma formation is not ruled out. 

Next we discuss the results obtained with in-medium meson masses. The separate contributions from the decay 
of mesons and baryons, and two-body processes are shown in Fig. 9. There is not much change in the meson and 
baryon decay contributions. With dropping meson masses, the abundances of p and ai mesons increase ||^,Q. This 
is, however, compensated by the decreases in their radiative decay widths (see Figs. 2 and 3). The comparison of 
total thermal photon spectra with the experimental data is given in Fig. 10. In the case of dropping meson masses, 
our results do not exceed the upper bound of the WA80 data. 

Finally, in Fig. II we show the ratio of total (sum of thermal and background) photon to the background photon. 
In both and free meson mass and in-medium meson mass cases, the thermal photon accounts for less than 5% of all 
the single photons, in agreement with the experimental observation, and the simple estimation by Tserruya in Ref. [^j. 
An important difference between dilepton and real photon measurements is that in the former case, the contamination 
from Tr° —> ye“'"e“ can, in principle, be removed since these dileptons are restricted to masses below while in the 
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real photon case, the —> 77 contributes to photons with any momentum. The subtraction of the background, and 

hence the measurement of interesting signals, are much more subtle in the real photon case j^. 


V. PHOTON PRODUCTION: COMPARISON WITH OTHER CALCULATIONS 

The WA80 single photon data have been looked at by several groups, mostly based on hydrodynamical models 
p2|- ^^5| - |3^ , ^ , ^ Here we compare our results that are based on the relativistic transport model with some of these 
results. In Fig. 12, we show our results together with those from and |Q. In Ref. ||^ different equations of 
state with and without quark gluon plasma formation were considered. Shown in the figure by dotted line is their 
results based on EOS H. This equation of state included about the same number of hadron degrees of freedom as in 
our transport model. It is very interesting to see that their results are very similar to ours, although the dynamical 
models used are quite different. In Ref. [pq two different equations of state, one with quark gluon plasma formation 
and one with pure pionic gas, were considered. Shown in the figure is their results based on the pionic gas equation 
of state. Because of very limited number of hadronic degrees of freedom in this equation of state, they needed a very 
high initial temperature to account for the final observed hadron abundances. This led to large photon yield. 

In Refs. Mm- Steele et al. studied photon as well as dilepton production at SPS energy based on the master 
equation formalism. They have considered two scenarios, one with pions only and the other includes nucleon as well. 
In order to make a fair comparison, we show in the upper window of Fig. 13 the photon spectra obtained in our 
calculation with the decays of p and ai and two-body processes of tttt ^ py and irp Try, and the results of Ref. p^ ] 
excluding nucleons. We can see that their results are very similar to ours. In the lower part of Fig. 13 we compare 
our results including baryon resonance decay with those of Ref. that included nucleon effects. The results are 
again quite consistent with ours. Note that the results of Ref. is very sensitive to the initial nucleon density. The 
results shown in the lower part of Fig. 13 were obtained with an initial nucleon density Pn = 0.7po, which is quite 
appropriate for S-|-Au collisions at SPS energies. 


VI. SUMMARY 

In summary, we studied single photon spectra in central S-|-Au collisions at SPS energies using the relativistic 
transport model that has been used to study dilepton spectra in the same reactions. We included photons from the 
background sources of 7 r° and p decays, as well as thermal sources such as meson decays, decays of baryon resonances, 
and two-body processes. We found that more than 95% of single photons come from the decays of 7 r° and 77 . The 
thermal photons account for only less than 5% of all single photons, in agreement with the experimental observation 
made by the WA80 and CERES collaborations. We compared our thermal photon spectra with the experimental 
upper bound extracted by the WA80 collaboration. It is seen that in both the free meson mass and in-medium meson 
mass cases, our results do not exceed the experimental upper bound. This indicates that our model can explain both 
the enhancement of low-mass dileptons and the lack of signal in the real photon measurement. Our results are also 
in agreement with the conclusions of hydrodynamical calculations of Refs. that included a sufficient number 

of hadron degrees of freedom. Therefore the WA80 single photon data do not necessary imply the phase transition to 
the quark gluon plasma, as claimed in Refs. based on a limited number of hadron degrees of freedom. 

Finally, we note that the elementary radiative decay widths and photon prodcution cross sections used in this study 
have been taken from different models. In principle, most of them can be evaluated consistently within the same 
model, such as the hidden gauge theory of |Q. In this model, one can also study the effects of the change of gauge 
coupling constant on photon production in hot and dense matter. Work in this direction is in progress and will be 
reported elsewhere ^ . 


We thank M.A. Halasz, C.M. Ko, J.V. Steele and I. Zahed for discussions, and I. Tserruya for a careful reading of 
the manuscript and useful suggestions. This work is supported in part by the Department of Energy under grant No. 
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Figure Captions 

Fig. 1: Pseudorapity distribution of charged hadrons (upper window) and pion transverse momentum spectra 
(lower window) in central S+Au collisions at 200 AGeV. 

Fig. 2: Decay width for p —> TTTry as a function of rho meson mass. The open circle is experimental data for . 

Fig. 3: Radiative decay widths of uj (upper window) and oi (lower window) mesons as a function of their masses. 
The open circles are experimental data. 

, for three different masses. 

, for three different masses. 

Fig. 6: Backgroud single photon spectra from and ry decays in central S+Au collisions. 

Fig. 7: thermal single photon spectra from p, uj, rj' and ai decays (left window), from N*, A* and decays 
(middle window), and from tttt and np scattering (right window) in central S+Au collisions. 

Fig. 8: thermal single photon spectra meson decay (dotted line), baryon decay (short-dashed), and two-body scat¬ 
tering (long-dashed) in central S+Au collisions. The solid circles give the upper bound from the WA80 collaboration 


Fig. 4: Cross sections for tttt ^ py from Ref. 
Fig. 5: Cross sections for irp Try from Ref. 


\M 


Fig. 9: Same as Fig 7, the results with in-medium meson masses. 

Fig. 10: Same as Fig. 8, the results with in-medium meson masses. 

Fig. 11: The ratio of total photon spectra to the backgroud photon spectra in central S+Au collisions. The solid 
and dashed histograms are the results with in-medium and free meson masses, respectively. The solid circles are the 
experimental data from the WA80 collaboration . 

Fig. 12: Comparison of our results with those of Ref. [p4| ( dotted line) and Ref. [Q (long-dashed line). 

Fig. 13: Comparison of our results with those of Ref. ||22|] . Upper window: pions only, lower window: pions and 
baryons. 
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